home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlacon.z / dlacon
Encoding:
Text File  |  2002-10-03  |  3.5 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAACCCCOOOONNNN((((3333SSSS))))                                                          DDDDLLLLAAAACCCCOOOONNNN((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLACON - estimate the 1-norm of a square, real matrix A
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE DLACON( N, V, X, ISGN, EST, KASE )
  13.  
  14.          INTEGER        KASE, N
  15.  
  16.          DOUBLE         PRECISION EST
  17.  
  18.          INTEGER        ISGN( * )
  19.  
  20.          DOUBLE         PRECISION V( * ), X( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      DLACON estimates the 1-norm of a square, real matrix A. Reverse
  37.      communication is used for evaluating matrix-vector products.
  38.  
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      N      (input) INTEGER
  42.             The order of the matrix.  N >= 1.
  43.  
  44.      V      (workspace) DOUBLE PRECISION array, dimension (N)
  45.             On the final return, V = A*W,  where  EST = norm(V)/norm(W) (W is
  46.             not returned).
  47.  
  48.      X      (input/output) DOUBLE PRECISION array, dimension (N)
  49.             On an intermediate return, X should be overwritten by A * X,   if
  50.             KASE=1, A' * X,  if KASE=2, and DLACON must be re-called with all
  51.             the other parameters unchanged.
  52.  
  53.      ISGN   (workspace) INTEGER array, dimension (N)
  54.  
  55.      EST    (output) DOUBLE PRECISION
  56.             An estimate (a lower bound) for norm(A).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAACCCCOOOONNNN((((3333SSSS))))                                                          DDDDLLLLAAAACCCCOOOONNNN((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      KASE   (input/output) INTEGER
  75.             On the initial call to DLACON, KASE should be 0.  On an
  76.             intermediate return, KASE will be 1 or 2, indicating whether X
  77.             should be overwritten by A * X  or A' * X.  On the final return
  78.             from DLACON, KASE will again be 0.
  79.  
  80. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  81.      Contributed by Nick Higham, University of Manchester.
  82.      Originally named SONEST, dated March 16, 1988.
  83.  
  84.      Reference: N.J. Higham, "FORTRAN codes for estimating the one-norm of a
  85.      real or complex matrix, with applications to condition estimation", ACM
  86.      Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
  87.  
  88.  
  89. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  90.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  91.  
  92.      This man page is available only online.
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.